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Considering model Hamiltonians that respect the symmetry properties of the pnictides, it is 
argued that pairing interactions that couple electrons at different orbitals with an orbital-dependent 
pairing strength inevitably lead to interband pairing matrix elements, at least in some regions of the 
Brillouin zone. Such interband pairing has not been considered of relevance in multiorbital systems 
in previous investigations. It is also observed that if, instead, a purely intraband pairing interaction 
is postulated, this requires that the pairing operator has the form A^(k) = /(k) a ^cft_ k a ^ 

where a labels the orbitals considered in the model and /(k) arises from the spatial location of 
the coupled electrons or holes. This means that the gaps at two different Fermi surfaces involving 
momenta k_F and k^ can only differ by the ratio /(ki?)//(k^) and that electrons in different orbitals 
must be subject to the same pairing attraction, thus requiring fine tuning. These results suggest that 
previously neglected interband pairing tendencies could actually be of relevance in a microscopic 
description of the pairing mechanism in the pnictides. 

PACS numbers: 75.50.Pp, 71.10.Fd, 72.25.Dc 



I. INTRODUCTION 

The discovery of superconductivity in the iron-based 
pnictidesi~— has opened a new active direction of re- 
search in the quest to understand high critical tempera- 
ture (high-T c ) superconductors. Experiments arc show- 
ing that the pnictides share several properties with the 
high-T c cuprates, such as the order of magnitude of the 
critical temperature) 1 - the existence of magnetic order in 
some of the parent compounds^— and a possible ex- 
otic pairing mechanism^ However, there arc differences 
in several aspects as well: the parent compound is a 
(bad) metal instead of a Mott insulator^— and sev- 
eral orbitals, as opposed to only one, have to be con- 
sidered in order to reproduce the Fermi surface, which 
consists of hole and electron pocketsJ^— In addition, 
while clear experimental evidence and theoretical calcu- 
lations indicate that the pairing state in the cuprates is 
nodal and has d-wave symmetry^ the properties of the 
pairing operator in the pnictides have not yet been estab- 
lished. Experimentally, several angle resolved photocmis- 
sion (ARPES) studies^— show constant nodcless gaps 
on all Fermi surfaces (FSs), but evidence for the exis- 
tence of nodal gaps has been reported in many transport 
measurements as welli^— It has been argued that the 
properties of the gap may be material dependent or that 
a nodal gap may be rendered nodeless by disorder^ 7 - but 
a consensus has not been reached. 

The goal of this paper is to understand the constraints 
that symmetries and the number of active degrees of free- 
dom in the pnictides impose on the possible pairing op- 
erators. We will consider a five orbital model that re- 
tains the five d orbitals of each of the two Fe atoms in 



the unit cell of the FeAs planes^ Employing mean- field 
approximations we will discuss comparisons with results 
obtained in models with thread and twc42r— orbitals 
that can be studied numerically. 

The organization of the paper is as follows: in Sec- 
tion [IT] we present the five-orbital Hamiltonian written in 
terms of SU(5) 5x5 matrices to properly identify its sym- 
metries. The pairing operators are explicitly discussed in 
Section IIIII Section [TV] is devoted to a discussion of the 
results, while conclusions are provided in Section IVl 



II. MODEL 

To construct the possible pairing operators allowed by 
the lattice and orbital symmetries, we will follow the ap- 
proach already described in detail for the case of a three- 
orbital models The first step involves the rewriting of 
the tight-binding Hamiltonian for the five-orbital model, 
for instance as provided in Ref. [38l in terms of some of 
the 25 5 x 5 matrix generators of C/(5) D U(l) x SU(5), 
which consist of the unit 5x5 matrix and 24 5 x 5 matrix 
generators of S£/(4). Since most of the elements of these 
matrices are zero, instead of providing for the reader this 
large number of matrices explicitly, it is more convenient 
to simply describe them verbally. For i = 1 up to 8, the 
Xi matrices are given by the well-known eight Gcll-Mann 
matrices^ in the upper left-hand corner while the rest of 
the elements are equal to zero. The matrices with i = 9 
(1,4), 11 (1,5), 13 (2,4), 15 (2,5), 17 (3,4), 19 (3,5), 
and 21 (4, 5) are symmetric and have only two elements 
equal to 1 while the rest of the elements are 0, with 
indicating the position of one of the non-zero elements 
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(by symmetry the other element can be identified). The 
matrices with i = 10 (1,4), 12 (1,5), 14 (2,4), 16 (2,5), 
18 (3,4), 20 (3,5), and 22 (4,5) are Hermitian and have 
two elements equal to ±i and the rest are 0, with (i, j) in- 
dicating the position of the non-zero element equal to —i. 
Finally, A23 is diagonal with nonzero elements (4,4) and 
(5,5) equal to 1 and -1, and A24 is diagonal with nonzero 
elements (1,1) = (2,2) = (3,3) = (4,4) = and 

The Hamiltonian then becomes 

ff TB (k) = £4 )(T C k $ k)<T , (l) 

k,<r 

where $ t a = (4(k), 4(k), 4(k), 4(k), 4(k)) ff and 

24 

C k = ^a 4 A l , (2) 

where the coefficients are presented in Table |U The 
functions appearing in this Table are provided in the 
Appendix. The on-site energies for each orbital are given 
by ei = e 2 = 0.13 cV, e 3 = -0.22 cV, e 4 = 0.3 eV, and 
£5 = —0.211 eV.— In addition, = + e^, and the 
chemical potential is 0. The index-to-orbital correspon- 
dence is the following: (1) xz, (2) yz, (3) x 2 — y 2 , (4) xy, 
and (5) 3z 2 -r 2 ^ 

The symmetry operations that leave invariant the Fe- 
As planes can be mapped on the elements of the point 
group -D4 H 43 ' 44 and, thus, the Hamiltonian has to re- 
main invariant under all the operations of this group, 
which means that it must transform according to the ir- 
reducible representation A ig . This allows us to assign an 
irreducible representation to each of the A 2 ; matrices and, 
thus, classify possible pairing operators according to their 
symmetr y 39 i It is important to realize that due to the 
strong hybridization among the orbitals, evident in all 
the terms in Eq. [T] where £y with i 7^ j appear, a proper 
characterization of the pairing operators by symmetry 
cannot be accomplished using only the band representa- 
tion. 



III. PAIRING OPERATORS 

In order to construct the allowed pairing operators in 
multiorbital models the symmetries of the spatial, spin, 
and orbital contributions needs to be considered. Ignor- 
ing the orbital symmetry may lead to problems similar to 
those encountered in the early days of the study of mag- 
netism when the spin contribution to the electronic wave 
functions was not included. Note that from the point of 
view of the orbital symmetry some of the previous ef- 
forts on superconductivity in multiorbital systems, such 
as in Ref. l46l are effectively dealing with non-hybridized 
"s-orbitals" . Thus, those results cannot be straightfor- 
wardly applied to a system with strongly hybridized non- 
s orbitals. In fact, it was observed that the general form 



TABLE I: Coefficients for the Ai matrices in Eq.[2]and the ir- 
reducible representation of according to which they trans- 
form. The £ij are provided in the Appendix. 
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of a spin-singlet pairing operator in this case is given by 4 ^ 

At(k) = / i (k)(A i )^(4, a ,/- k ^,,-4,^/-k,«J> (3) 

where a sum over repeated indices a and /3 is implied; the 
operators d have been defined above, and /(k) is the form 
factor that transforms according to one of the irreducible 
representations of the crystal's symmetry group^ The 
form factors that will be considered in this work and their 
corresponding irreducible representations of (accord- 
ing to which they transform) are presented in Table [TTJ 
The index i in Eq. ([3]) indicates that different form factors 
may be needed if matrices Aj with different symmetries 
are combined. 

As in the case of the twc4^ and three 3 ^ orbital models, 
the properties of the pairing operators can be studied un- 
der a mean-field approximation. We need to remember 
that the five-orbital model is defined in terms of a pseu- 
docrystal momentum k in an extended Brillouin zone^ 
In terms of the real momentum, the unit cell of the Fe-As 
planes contains two Fe ions and, thus, the band struc- 
ture is composed of 10 bands in the reduced, or folded, 
Brillouin zone (BZ). This can be observed in Fig. Ufa) 
where the spectral function A(k, ui) is shown for the non- 
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TABLE II: Form factors /(k) for pairs up to distance (1,1) 
classified according to their symmetry under Da^ operations. 
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interacting case along high symmetry directions in the 
folded BZ where the 10 bands are clearly seen. Thus, for 
each real momentum q there are 10 bands. Five of them 
are the bands of the five-orbital model for the pseudocrys- 
tal momentum k = q and the other five are obtained from 
the same model by setting the pseudocrystal momentum 
to k = q + Q, where Q = (it, it). Since in terms of the 
real momentum, the basis of the five-orbital model is ex- 
panded by states with momentum q for orbitals xz and 
yz and momentum q + Q for the other three orbitals, 
this fact needs to be taken into account when pairing op- 
erators are being constructed. If only intraorbital pairing 
operators are considered it is sufficient to build a 10 x 10 
Bogoliuvov-deGennes (BdG) matrix, but for interorbital 
pairing between electrons in orbitals 1, 2 with electrons 
in orbitals 3, 4, or 5, it is necessary to construct a 20 x 20 
BdG matrix since the 10 x 10 matrix only allows to con- 
sider intraorbital pairs with pseudocrystal momentum Q 
rather than 0.— 

The 10 x 10 BdG Hamiltonian is given by: 

H BdG = j2*K F ^ (4) 

k 

with the definitions 

vt^ — (eft d^ (ft d^ 

— V u k,l,t' "k,2,f U k,3,f "k,4,f "k,5,f 
d-k,l,l, ^-k,2,4,3 ^-k,3,4,^-k,4,4,, eLk,5,4,)j (5) 

and 

where each element represents a 5 x 5 block with iZrB(k) 
given by Eq. (JTJ) and 

P(k) ai/3 = Vf(\n)(\i) a ,p, (7) 

where V = VqA is determined by the product of an un- 
known pairing strength Vo and a parameter A that arises 
from minimizing the mean-field equations, as already ex- 
plained in previous literature4£ 
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FIG. 1: (Color online) The intensity of the points represents 
the values of the spectral function A(k,cj) for the five-orbital 
model with pairing interaction (a) V = 0; (b) V = 0.2, and 
for the 5± pairing operator given in the text. The results are 
shown in the reduced BZ. 



A. Intraorbital Pairing Operators 

1. Purely Intraband Pairing Operators 

The first issue we want to address is what is the form 
of the pairing operators resulting from purely intraband 
pairing interactions. The motivation is given by the fact 
that in standard BCS theory, the pairing occurs between 
particles with momentum equal in magnitude but oppo- 
site direction at a common Fermi surface. For this reason, 
even in multiorbital systems, it has been expected that 
pairing attraction should involve particles in the same 
band4£ To determine whether a pairing operator consists 
of purely intraband matrix elements, we need to trans- 
form Eq.[B]from the orbital to the band representation in 
which Htb is diagonal. The transformation is given by 
#Band(k) = W(k)H TB {k)U{k), where U(k) is the uni- 
tary change of basis matrix and (k) is the transpose 
conjugate of U(k). Since U is unitary, for each value 
of k, J2i(Ui,j)*Ui, m = Y,i{Uj,i)*U m ,i = <W Then, 
ff^p = G'HytpG where G is the 10 x 10 unitary ma- 
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trix composed of two 5x5 blocks given by U . Then, The most general form of a purely intraband pairing is 

given by 

rr/MF _ / -ffBand(k) ffe(k) \ , . 

k "1 4(k) -ifBand(k) ) ' ^ 

with 

P B (k) = [/- 1 (k)P(k)[/(k). (9) 
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where £j(k) are the eigenvalues of iJxB(k), while Aj(k) 
denotes the band and momentum dependent pairing in- 
teractions. Notice that if A^ = Ao in Eq.0, -P(k) in 
Eq.® is proportional to the identity matrix and, thus, 
P B (k) is diagonal with Aj(k) = A(k) = Vf{k) for all 
j. This indicates that the intraorbital operator that 
pairs electrons in each orbital with equal strength, gives 
rise to a purely intraband pairing operator where elec- 
trons in each band are subject to an identical pairing 
attraction. In this case, diagonalizing Eq. ([10[) we ob- 
tain Ej(k) = ^(k) 2 + |A(k)| 2 . Then, at k F where 
ej(kp) = we obtain Ej(k F ) = |A(k^)|. This means 
that the superconducting gaps at the FS determined by 

different bands must satisfy that A'(k^)/A(ki?) = J7^f] 
where k^ and k^ represent the Fermi momentum of two 
FS defined by two different bands. In the case of the stan- 
dard low-temperature BCS pairing /(k) = 1, implying 
that momentum independent gaps of equal magnitude 
should open in all the FSs. However, for the pnictides 
it is believed that a non-BCS interaction provides the 
source of pairing*^ Then, we must consider the case in 
which /(k) 7^ 1. Notice that since Ao transforms accord- 
ing to Ai g , the symmetry of this purely intraband pairing 
operator is then totally determined by the symmetry of 
/(k). 

2. Single Gap: The 5± pairing operator 

In the particular case in which f(k) — cos k x cos k y , 
P(k) transforms according to Ai g and it represents the 



simplest form of the well known 5± pairing stated - — 
which will be characterized by gaps of approximately the 
same magnitude in the electron and hole pockets if there 
is good (it, 0) and (0, ir) nesting. Note that the gap on 
the two hole FSs would differ only if they have consid- 
erably different Fermi momentum which is not the case 
in the pnictides since both hole pockets are very close 
to each other. In the five-orbital model we found that 
r = jft4 ~ 1 for the external hole and electron pockets, 
while r ~ 0.9 for the external and internal hole pockets. 
This means that the gaps in the two hole pockets would 
have a difference of only about 10%. The two Fermi sur- 
faces would have to be much more separated from each 
other in order to develop the experimentally observed 
difference of 50% reported via ARPES experiments.— 

In Fig. [2] we show the gap as a function of the angular 
position (from the x to the y axis) along the four Fermi 
surfaces (two electron and two hole pockets) for the five- 
orbital model in the folded BZ. Note that the values of the 
gap beyond $ = 7r/4 can be obtained by symmetry from 
the values shown in Fig. [2J The interaction V = 0.02 
was chosen in order to obtain a quantitative match with 
experimentally reported values of the gaps that range be- 
tween 1 to 20 meVi2£ As expected, it can be seen that the 
gaps at the external hole and electron pockets (dashed, 
continuous, and dotted lines) almost overlap with each 
other and have a magnitude of about 7 meV while at 
the internal hole pocket (dotted-dashed line) the gap is 
around 8 meV, i.e. only about 10% larger. Note that 
the momentum dependence observed here likely will be 
negligible when considering experimental uncertainties. 
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Thus, these results would be in good agreement with the 
ARPES measurements reported in Ref. [U if we assume 
that the 50% smaller gap arises in a third hole pocket, 
rich in x 2 — y 2 orbital content, that is not present in this 
five-orbital model. 

The mean-field calculated spectral functions are shown 
in Fig. QJb) along high symmetry directions in the re- 
duced BZ. In the figure, 1/—0.2 is used 5 ^ and it can 
be seen that the gaps, that are very similar in magni- 
tude, have opened on all the original FSs. As expected, 
a numerical inspection of the whole BZ shows that there 
are no nodes. It can be observed that band distortions 
and "shadow" (or Bogoliubov) band spectral weight^ 5 - 
develop only in a small region around the chemical po- 
tential. 

If a different form factor such as /(k) = cos k x — cos k y 
is considered, the pairing operator would still be purely 
intraband but with symmetry -Big. While now nodes 
would appear on the FS because the pairing interaction 
vanishes for k x = ±k y , the gaps on the different FSs still 
only will differ by the ratio of the form factors at the 
respective Fermi momenta. 
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become non-diagonal in the band representation, i.e., 
Af = {7 _1 (k)Aji7(k) is non-diagonal at least for some 
values of k in the BZ. This is true even for A^ matrices 
that are diagonal in the orbital representation such as 
A3, As, A23, and A24. Thus, all these purely intraorbital 
pairing operators lead to interband matrix elements in 
their band representation. Note that to generate an 
intraorbital pairing operator that couples electrons in 
different orbitals with arbitrary strengths we need to 
consider linear combinations of these diagonal pairing 
matrices using form factors that provide a well defined 
symmetry, since we assumed, guided by numerical simu- 
lations in models with two and three orbitals, that the 
pairing operator connects non-degenerate ground states 
that transform according to one of the one-dimensional 
irreducible representations of D^. 

We will consider a linear combination of the four in- 
traorbital pairing operators allowed by symmetry that do 
not require the pairing to be the same in all orbitals. Ao 
is excluded because we already know that it does not pro- 
duce interband pairing and, thus, it will not contribute 
to off-diagonal elements in the band representation of the 
pairing operator. The linear combination to be consid- 
ered is given by: 

A = aiA 8 /i(k) + a 2 A 2 3/2(k) + 0^24/3 (k) + a 4 A 3 /4(k), 

(11) 

where /i(k) are form factors chosen such that all the 
terms in A transform according to the same irreducible 
representation. Notice that Ag, A23, and A24 transform 
according to Ai g , but A3 transforms according to Bi g . 
Thus, for example, a nearest-neighbor pairing operator 
that transforms according to Ai s will be obtained by 
setting /i(k) = / 2 (k) = / 3 (k) = cosfc : 
/ 4 (k) = cost 
becomes 



cos k y and 
cos k y . In the band representation A 



FIG. 2: (Color online) The gap at the internal-hole (dashed- 
dotted line), external- hole (dashed line), external-electron 
(continuous line) and internal-electron (dotted line) Fermi 
surfaces for the S± pairing operator with V = 0.02. Results 
are shown as a function of the angle <E> between and 7r/4 mea- 
sured with respect to the k x axis in counterclockwise [clock- 
wise] direction for the hole [electron (at X)]-pockets. The 
results are in the folded Brillouin zone. 



3. Multiple Gaps 

Since several experimenta l 21 ' 23 ' 33 ' 56 ' 57 and 
theoretica l 58 ' 59 efforts have reported the existence 
of at least two independent gaps at the Fermi surfaces 
of some pnictides, the next issue to consider is whether 
there is any other pairing state allowed by symmetry that 
is purely intraband and able to generate independent 
gaps in at least one of the FSs. It can be shown that, 
with the exception of Ao, all the other 24 Xi matrices 



A 6 = C/- 1 (k)Ai7(k). 



(12) 



If the pairing is purely intraband then Af, has to be diag- 
onal. But from the orthogonality properties of U we see 
that a non-diagonal element (Ab)ij can only be zero if 
A^j is the same for all i. However, this is only true for the 
pairing operator proportional to Ao and we have proven 
numerically that only for otj = in Eq.[TT]is A(,(m, n) = 
for m =/= n. Thus, intraorbital pairing with orbital de- 
pendent strength leads to interband components of the 
pairing interaction^ 

Then, we conclude that if the pairing interaction is 
purely intraband the pairing mechanism should be as- 
sociated to a degree of freedom that couples to the five 
different orbitals with the same strength. This means 
that the gaps in all different FSs must be related, i.e., 
the symmetry does not allow unrelated gaps in this case. 
Conversely, if independent gaps are observed on different 
FSs the symmetry of the highly hybridized orbitals in- 
dicates that interband interactions would be present at 
least in some regions of momentum space, as it will be 
discussed below ^ 
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FIG. 3: Band structure of the five-orbital model along high 
symmetry directions in the reduced Brillouin zone. The irre- 
ducible representations of Z?4h(r) and D2h(X), characterizing 
the bands at high symmetry points in the reduced BZ, are in- 
dicated. For more details see also the discussion in Sec. IIVI 



Some ARPES data3i indicate a superconducting state 
with momentum independent gaps with value A at the 
hole and electron pockets and A/2 at a third hole pocket 
not present in the five-orbital model, at least with the pa- 
rameters used here and in Ref . l38l. However, a slight mod- 
ification of the parameters, without affecting the sym- 
metry of the Hamiltonian,— would create an additional 
hole pocket around r of x 2 — y 2 character. This is the 
state with symmetry at the T point shown in Fig. [3] 
According to this figure, close to the hole FSs around 
r the pairing matrix should be diagonal with elements 
A (A/2) for the bands labeled E g (B| g ) at T, which in 
the unfolded BZ corresponds to (i) a diagonal pairing 
matrix at k « k^. [h denoting the hole Fermi surface) 
with at least two diagonal elements equal to A, and (ii) 
a diagonal pairing matrix at k « k^ + Q with at least 
one diagonal element equal to A/2. As just discussed, 
these two different values of a momentum-independent 
gap cannot arise from purely intraband pairing interac- 
tions in a highly hybridized system. 



4- Nodeless Gaps with Interband Pairing Matrix Elements 

Since several experimental ARPES studies of the pnic- 
tides appear to indicate that nodeless gaps open at all the 
FSs in the superconducting state^Sr— we will now iden- 
tify the pairing operators that produce nodeless gaps that 
are allowed by the symmetries of the five-orbital model. 

In a previous study of a three-orbital model we found 
that, in addition to the Si pairing operator, there was 
another nodeless pairing operator with both intra and 
interband matrix elements that was called Sib-— It leads 
to intraorbital pairing with different strength for the or- 
bitals xz/yz and xy, and it transforms according to Ai g . 



In the five-orbital model a similar result has been ob- 
tained. In fact, we have found several linear combinations 
with Ai g symmetry of the intraorbital pairing matrices 
Ai with i = 0, 8, 23, and 24 that provide nodeless gaps. 
In these pairing states, /(k) = 1 or cosk x cosk y and the 
pairing interaction is not the same in the different or- 
bitals. The pairing matrix has the form 



/ Vl \ 

Vi 

V 2 

v 3 o 

V Vi J 



(13) 



where Vi denote the different pairing strengths. Exam- 
ples of parameter sets for which nodeless gaps are found 
are {Vx,V u V 2 ,V z ,V A ) = V(l, 1, 0, 2, 0), 7(1,1,0,1,0), 
and V(l, 1, 0, 0, 0). These operators pair electrons in the 
orbitals that contribute the most to the FS, but it is im- 
portant to notice that they are not diagonal in the band 
representation. 

The momentum dependence of the gaps at the FS 
is shown in Fig. |4(a)| for the Sib pairing operator 
with /(k) = cosk x cosk y and (Vi, V\, V 2 , V3, V4) = 
(0.01, 0.01, 0.0, 0.02, 0.0)/V6, as special case. To repro- 
duce experimental values for the gap2£ V = 0.01/-\/6 
has been chosen. It can be observed that in this case 
the gaps for the internal hole and electron pockets in 
panel (a) (dashed-dotted and dotted lines) have a neg- 
ligible momentum dependence and an average value of 
about 5 meV. The external hole pocket has a smaller gap, 
about 3.5 meV, with a very weak momentum dependence 
(dashed line). The strongest dependence with momen- 
tum occurs in the external electron pocket (continuous 
line) with the gap reaching a value of about 7.5 meV at 
the point along the x-axis where the electron pocket is 
the closest to the hole pockets and has pure xy orbital 
character. As the point where the two electron-like FSs 
intersect each other (4> = 7r/4) is approached, the ex- 
ternal electron-pocket acquires xz/yz character and the 
magnitude of the gap converges to the almost momen- 
tum independent value of 5 meV that characterizes the 
internal electron-pocket. Thus, this is a case in which a 
direct measurement of the gaps would indicate the pres- 
ence of three "independent" gaps. The gaps in the in- 
ternal electron and hole pockets behave as expected in 
the S± pairing scenario ^2r— while the external hole and 
electron pockets show a different gap value. While this 
would agree with experimental results that favor mul- 
tiple gap s 21 ' 23 ' 33 ' 56 " — it is important to realize that this 
pairing operator contains interband matrix elements that 
must be incorporated in the description of a possible pair- 
ing mechanism. 

Modifying Vi we can tune the relative gap values, ob- 
taining three almost momentum independent gaps with 
values 3.5 meV, 5.5 meV, and 7.5 meV (in the external- 
hole, internal-hole, and internal-electron pockets, respec- 
tively) as shown in Fig. |4(b)| In addition, a gap with 
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FIG. 4: (Color online) The gap at the internal-hole (dashed- 
dotted line), external- hole (dashed line), external-electron 
(continuous line), and internal-electron (dotted line) Fermi 
surfaces for the nodeless pairing operator Sib with (a) 
(V u V 1 ,V 2 ,V 3 ,V i ) = (0.01, 0.01, 0.0, 0.02, 0.0)/V6 and (b) 
(Vi, Vi, V 2 , V 3 ,Va) = (0.005,0.005,0.0,0.02,0.0). Results are 
shown as a function of the angle $ between and 7r/4 mea- 
sured with respect to the k x axis in counterclockwise (clock- 
wise) direction for the hole [electron (at X)]-pockets. The 
notation is as in Fig. [2] 



ets that intersect at two points) then the formation of 
interband pairs, as it was described in Ref. H3, could 
occur. In addition, numerical Lanczos studies of a two- 
orbital model for the pnictides suggest that an intcror- 
bital pairing state with i?2 g symmetry, with interband 
components, could be the favored pairing state in the in- 
termediate Hubbard U rcgime42r— For this reason, here 
results will be presented for some intcrorbital pairing op- 
erators that have interband attraction, at least in some 
regions of the BZ. The effects of this interaction will be 
addressed later in Sec. HVl 
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FIG. 5: (Color online) The gap at the internal-hole (dashed- 
dotted line), external-hole (dashed line), external-electron 
(continuous line), and internal-electron (dotted line) Fermi 
surfaces for the interorbital pairing operator B 2g with V = 
0.012/v / 5- Results are shown as a function of the angle $ 
between and n/4 measured with respect to the k x axis in 
counterclockwise [clockwise] direction for the hole [electron 
(at X)]-pockets. The notation is as in Fig. [2] 



First let us consider the B2 g pairing operator men 
tioned in the previous paragraph. It is given by 



momentum dependence ranging between 15 meV and 7.5 
meV appears at the external electron pocket which has 
the largest xy contribution. Thus, this interaction would 
lead to four apparently independent nodeless gaps. 

B. Interorbital Pairing Operators 

As discussed above, due to the strong hybridization of 
all five d orbitals we have verified that all the intraor- 
bital pairing operators allowed by the lattice and orbital 
symmetries of the pnictides lead to interband pairing in- 
teractions if Xi ^= Ao- Interband pairing has always been 
considered unlikely in BCS theory 4 ^ because the pairing 
attraction acts in a very narrow energy range around the 
FS. However, in a system in which two FSs formed by 
different bands are very close to each other (e.g. the two 
hole pockets in the pnictides that cannot be distinguished 
in ARPES experiments ) 21 ' 23 ' 24 or the two electron pock- 



2 

A k g ( k ) =V(cosk x + cos k y ) 4,a,t d -k,/u- ( 14 ) 

For this operator, the structure of the gap strongly de- 
pends on the value of the pairing attraction V. For values 
of V compatible with the order of magnitude of the gaps 
reported from experiments in the pnictides, which are of 
the order of meV, the behavior of the gap at the different 
FSs is shown in Fig. [5] The gap on the two hole pockets 
presents nodes along the x (3? = 0) and y = ir/2) 
axes (the latter not shown explicitly since the results are 
symmetric and thus A($) = A(7r/2 — $)). Both gaps 
are maximized along the diagonal direction ($ = 7r/4). 
Thus, it would resemble "e^-wave" behavior in experi- 
ments where only the node location but not the phase of 
the gap can be measured. The reason for the existence 
of the nodes is that along the x and y axes the pair- 
ing interaction is purely interorbital, but V is not strong 
enough to open a gap (a finite value of V is needed to 
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open a gap with purely intcrorbital pairing, as discussed 
in previous literature 4 ^) and that effect creates the node. 
By contrast, in Fig. EJa) we show the mean-field calcu- 
lated spectral functions along some high symmetry di- 
rections of the reduced BZ for the same pairing operator 
but with a larger V = 0.2. In this case, the gaps have 
opened on the hole-pocket FSs around the T point. But 
for this value of V the gap is of the order of 100 meV, 
i.e. too large compared with the experimental results for 
the pnictides. 

Another characteristic of the B 2g pairing operator is 
that the gap in the electron pockets presents nodes along 
the x and y axes (for the same reason than the hole pock- 
ets at weak V) but also at the points where the two elec- 
tron pockets cross, <I> = 7r/4, because /(k) = there. For 
V < 0.01 the gap is much smaller than the one on the 
hole-pockets, as shown in Fig. [5j As V increases, a gap 
opens along the x and y axes for the internal electron 
pocket, as it can be observed in Fig. [SJa) for V = 0.2. 
However, in this case nodes remain along these axes for 
the external electron pocket. The reason is that at these 
points the FS arises from bands that have mostly char- 
acter xz/xy or yz/xy^- and, thus, the operator Ag 2g , 
which couples electrons in the xz/yz orbitals, is not ef- 
fective at opening a gap along the momentum axis where 
one of the FSs has purely xy character. While this gap 
structure for the B 2g pairing state disagrees with ARPES 
measurements it is important to keep in mind that the 
existence of nodes in some pnictides has been reported 
by several groups using other experimental techniques. 
Thus, it is still possible that the surface that is actually 
tested by ARPES does not present the same behavior as 
the bulk in the pnictides. 

A natural generalization of the operator A^ 2 to in- 
clude the xy orbital was provided in Ref. HH. The ex- 
tended operator is a linear combination of (Eq.ll4p 
and 

sinfc a ^2 4,a,t d -k,/3j- ( 15 ) 

The total pairing operator has symmetry B 2g . Note that 
now the 10 x 10 matrix given in Eq. [4] will provide pairs 
with pseudocrystal momentum for the particles in or- 
bitals 1 and 2, but Q for one of the particles in orbital 
4 and the other in 1 or 2. Thus, we must consider the 
20 x 20 BdG matrix mentioned earlier and we have stud- 
ied the extended B 2g pairing with the two possible pseu- 
docrystal momenta and Q. This generalized form re- 
moves the nodes at the crossing point of the two electron 
pockets, since Ay is finite at those points. 

In Fig. Hth) we show the spectral function for V = 
V' = 0.2 for pairs with pseudocrystal momentum. It 
can be seen that isolated nodes at the electron pockets re- 
main only at the points along the momentum axis where 
they have purely xy character. However, for smaller val- 
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FIG. 6: (Color online) The intensity of the points represent 
the values of the spectral function A(k,cj) for the five-orbital 
model with pairing interaction V = 0.2, for the pairing oper- 
ator (a) B2 g and (b) extended f?2 g discussed in the text. The 
results are shown in the reduced BZ. 



ues of V = V that produce gaps in the meV range, we 
have found that nodes along the axis remain and the gaps 
resemble very much those obtained with the original B 2g 
pairing operator shown in Fig. [5j 



IV. DISCUSSION 

In theoretical studies of multiorbital systems it 
has been "traditional" to consider non-hybridized 
orbitals/bands4£ While this approach works for simple 
cases, in the pnictides the hybridization of the orbitals is 
strong. 14 ' 15 As a result, the energies in the band represen- 
tation €j(k) have several accidental degeneracies within 
the Brillouin Zone which means that the bands cross at 
several points and are, thus, very entangled^ 5 - This band 
entanglement is apparent in Fig. [31 where the eigenval- 
ues of £/tb are shown along high symmetry directions in 
the folded BZ. The irreducible representations character- 
izing the bands at some high symmetry points (such as 
r with symmetry and X with symmetry D 2 h) are 
indicated, and the star labels bands determined by states 
with pseudocrystal momentum k + Q. For example, the 
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band labeled Bi g at T, crosses two other bands along 
the r — X direction. Thus, the approach used in early 
studies by Suhl et ai, where a gap Aj is associated to a 
band with energy £j(k) in all the BZ, becomes ambigu- 
ous in this case. In a system with strongly hybridized 
orbitals, it is more reasonable to define pairing operators 
in the orbital representation because the orbital basis is 
globally well-defined at all points in the BZ, while the 
band assignation is local. In addition, as it was shown in 
the previous Section, the orbital representation provides 
the natural framework to construct the pairing operators 
that are allowed by symmetry. All the active degrees of 
freedom need to be incorporated in the pairing operators 
since otherwise the results can be misleading, as in the 
early studies of magnetization when the contribution of 
the spin to the spatial wave function was disregarded. 
Then, when intraband pairing operators are constructed 
their orbital content should be provided. 

Another important point is to understand the conse- 
quences of interband matrix elements arising from the 
symmetry of the pairing operators. As shown in Ref. \at\ , 
at points in the BZ where there are intra and interband 
elements in ffe(k) the pairs will be formed by electrons 
in the same band if the pairing interaction is weak or in- 
termediate, which seems to be the case in the pnictides. 
It is only at the very few points on the FS where PbO*-) 
has finite (zero) non-diagonal (diagonal) elements that 
the pairing attraction will be purely interband. At these 
points, if the pairing attraction is weak we would expect 
to observe nodes 4^ It is only when the two FS are very 
close (as for the two almost degenerate hole pockets), 
and when the pairing interaction is strong enough, that 
interband pairs would be possible^ 

As pointed out in the Introduction, the properties of 
the pairing operator in the pnictides have not been estab- 
lished yet and data obtained with different experimental 
techniques are in disagreement. ARPES results are in- 
terpreted as indicative of nodeless gaps with at most a 
very weak momentum dependence. Some groups have 
identified gaps with the same magnitude in the electron 
and one of the hole pockets,— ~— while others have found 
them to be different^ In addition, the gap in a third hole 
pocket, not of xz/yz character, is found to be different 
from the gap in the electron pockets and, thus, the exis- 
tence of two and, even three, gaps in these materials has 
been proposed. The symmetry arguments presented here 
indicate that truly momentum independent gaps would 
have to be equal (including the sign of the order parame- 
ter) in all Fermi surfaces. If a momentum dependence is 
allowed, then a pairing state such as the Si could arise. 
In this case the gap in the hole and electron pockets re- 
lated by nesting should be very similar (although with 
opposite signs in the order parameters) and the gap in 
the additional hole pocket would differ only by the ratio 
between the form factors /(k) at the locations of the two 
different hole-pockets. Thus, there should be only one 
gap in the sense that the coupling between the electrons 
and the interaction causing the pairing will not be orbital 



dependent and the band dependence is only due to the 
different Fermi momenta. 

We have also shown that a pairing interaction that 
couples with an orbital-dependent strength to the elec- 
trons leads to a pairing matrix that is not diagonal in 
the band representation. Some of these operators open 
weakly momentum dependent nodeless gaps in the dif- 
ferent FSs. While this kind of pairing operator could 
agree with some experimental results, it would be neces- 
sary to take into account the interband pairing interac- 
tions when developing the associated microscopic pairing 
mechanism. 

Note that the existence of nodes in the pnictides is in 
fact supported by transport experiments^^— in clear 
contradiction with ARPES. Thus, more experimental 
work is needed in order to clarify this issue. We have 
found a large variety of nodal pairing operators that re- 
spect the symmetry of the pnictides: this includes those 
proportional to the identity in the orbital sector, with 
their nodes arising from the zeroes in /(k), and also other 
pairing operators in which electrons in different orbitals 
are subject to different pairing strengths and, thus, give 
rise to interband terms in the pairing matrices. 

The Si pairing operator with symmetry A\ g appears 
to be the favored one in the literature; 4 ^— but numerical 
calculations in a two-orbital mode l 40 ' 41 while indicating 
that this pairing prevails in strong coupling, lead in the 
intermediate regime to a pairing state made of electrons 
in different orbitals and with symmetry #2g- Since the 
orbitals xz and yz are strongly hybridized forming the 
bands that produce the hole pockets, this kind of pair- 
ing could be possible. If realized, it would induce nodes 
on the hole and electron pockets and the pair formation 
would be much stronger in the hole than in the electron 
pockets. Interestingly, the gaps in the two hole pockets 
would be different (see Fig. [5]) which may be in agree- 
ment with some of the experimental results that indicate 
two nodal gaps^ 



V. CONCLUSIONS 

Summarizing, we have shown that in a model that re- 
tains the symmetry of the FeAs planes for the pnictides 
and considers the five d orbitals of the Fe ions, a purely 
intraband pairing operator can only result from an in- 
traorbital pairing interaction that affects electrons in the 
different orbitals with identical strength. In this case, the 
symmetry of the pairing operator is entirely determined 
by the form factor /(k) which depends only on the spa- 
tial location of the particles that form the pairs. As a 
result, gaps in different portions of the FS can differ only 
by the ratios of the form factors, i.e., two or more un- 
related gaps cannot occur. Conversely, multiple gaps as 
observed experimentally, or orbital dependent pairing at- 
tractions, would indicate interband pairing interactions 
at least in some regions of the BZ. Then, this feature 
should be incorporated in theoretical proposals for the 
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pairing mechanism. If there are special points in the BZ 
where the attraction is purely interband, then nodes or 
interband pairing will occur depending on whether the in- 
teraction is weak or strong. Experimental measurements 
of the gap magnitude indicate that the pairing attraction 
is weak, thus nodes rather than interband pairing would 
be expected. 

The present analysis suggests that if the pairing mech- 
anism is purely intraband, as assumed by many, then 
no unrelated gaps should occur in the different portions 
of the FS. Reciprocally, if it is experimentally confirmed 
the existence of two or more unrelated gaps this would 
point to the need to consider interband, in addition to in- 
traband, pairing interactions in any realistic microscopic 
description of the pairing mechanism. 

VI. ACKNOWLEDGMENTS 

This work was supported by the National Science 
Foundation grant DMR-0706020 and by the Division of 
Materials Sciences and Engineering, Office of Basic En- 
ergy Sciences, U.S. Department of Energy, M. D. ac- 
knowledges partial suppost from the DFG under the 
Emmy-Noether program. 

Appendix A: Expressions for the Tight Binding 
Hamiltonian 

£u = 2t lx cos k x + 2t l y cos k y + &b xy cos k x cos k y + 
2t xx (cos2k x — cos2k y ) + 4t xxy cos 2k x cosk y 
+4t xyy cos 2k y cos k x + 4:t xxyy cos 2k x cos 2k y (Al) 

£22 = %ty cos k x + 2t x cos k y + 4t xy cos k x cos k y — 
2t x \ (cos 2k x — cos 2k y ) + 4i** cos 2k x cos k y 
+4t xxy cos 2k y cos k x + 4t xxyy cos 2k x cos 2k y (A2) 

£33 = 2i 33 (cos k x + cos k y ) + 4t xy cos k x cos k y + 



2tf x (cos 2k x + cos 2k y ) (A3) 

£44 = 2i 4 , 4 (cos k x + cos k y ) + 4t X y cos k x cos k y + 

2t^ x (cos2k x + cos2k y ) 

+4i 44 y (cos 2k x cos k y + cos 2k y cos k x ) + 

4t xxyy cos 2k x cos 2k y (A4) 

£ 55 = 2tf (cos k x + cos ky ) + 2t x 5 x (cos 2k x + cos 2k y ) + 

4:t x 5 xy (cos 2k x cos k y + cos 2k y cos k x ) 

+4t x 5 xyy cos 2k x cos 2k y (A5) 

£12 = At x 2 y sin k x sin k y + 

At xxy (sin 2k x sin k y + sin 2k y sin k x ) + 

KlyySin2k x am2k y (A6) 

£13 = 2it x 3 sin k y + iit xy sin k y cos k x — 

iit xxy (sin 2k y cos k x — cos 2k x sin k y ) ( A7) 

£23 = 2it x 3 sin k x + 4it xy sin k x cos k y — 

4it xxy (sin 2k x cos k y — cos 2k y sin k x ) ( A8) 

£14 = 2it x 4 sin k x + 4it xy cos k y sin k x + 

iit x % y sin 2k x cos k y (A9) 
^24 = — 2H X 4 sin fc. y + iit xy cos fc x sin fc. y — 

kit x A xy sin 2fc y cos k x (A10) 
^15 = 2ii^ 5 sin k y — 4«i^ sin k y cos /c^ 

-Ait l xxyy sm2k y cos2k x (AH) 
£25 = —^H x b sin fcj; -I- 4it^ sin fc^ cos k y 

+4it x 5 xyy sm2k x cos2k y (A12) 
^34 = 4t^L(sin2fc s/ sinfc x — sin2fe x sinfcy) (A13) 
^35 = 2t x 5 (cosk x — cos k y ) + 

4t xxy (cos 2k x cos k y — cos 2k y cos fc^) (A14) 
^45 = 4t xy sin fej; sin k y + it xxyy sin 2^ sin 2k y (A15) 

The values of the hopping parameters are explicitly 
provided in Ref. HH. 
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